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O ■ Abstract 
O 
(N 

^ ! A formulation of quantum electrodynamics is given that applies to atoms in a 

I strong laser field by perturbation theory in a non-relativistic regime. Dipole 

' approximation is assumed. With the current wavelengths, squeezing can be 

proved to be negligible and then, just the linear term in the Hamiltonian 
. can be retained. The dual Dyson series, here discussed by referring it to the 

, Birkhoff theorem for singularly perturbed linear differential equations, can 

I be applied and a perturbation series obtained transforming the Hamiltonian 

' by a Pauli-Fierz transformation. But, if just few photons are present high- 

■ order harmonics cannot be generated. So, it is proven that odd high-order 

. harmonics only appear when the laser field is intense and one can substitute 

, the creation and annihilation operators by the square root of the mean num- 

ber of photons taken to be huge, the field retaining its coherency property as 
observed experimentally for harmonics. In this case, the Hamiltonian for per- 
ip5 ' turbation theory comes to the Kramers-Henneberger form. From this Hamil- 

tonian it is shown that just odd harmonics of the laser frequency contribute 



Q^. to the spectrum for a spherically symmetric potential. This contribution is 

dipolar when the free-electron quiver motion amplitude is larger than the 
atomic radius. For a Coulomb potential one has that the outer electron is 



' periodically kicked, and so a prove is given that the same should happen to 

^ Rydberg atoms in intense microwave fields. The distribution representing the 

kicking has a Fourier series with just odd terms. Using a modified Rayleigh- 
Schrodinger perturbation theory, it is shown that under the same condition 
of validity of the quiver motion amplitude to atomic radius ratio, the atomic 
wave function is only slightly modified by the laser field due to the way the 
energy levels rearrange themselves. This gives a prove of stabilization in the 
limit of laser frequency going to infinity. Then, perturbation theory can be 
applied when the product of the laser frequency and the square root of the 
ratio between the ionization energy and the ponderomotive energy, that is 
the Keldysh parameter, becomes smaller with respect to the shifted distance 
between the energy levels of the atom. 
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I. INTRODUCTION 



The availability of powerful sources of laser light has permitted in recent 
years the study of light-matter interaction in regimes where known approxi- 
mations fail. Indeed, a number of new phenomena has appeared as high-order 
harmonics of the laser frequency, ionization with a number of photons well 
above the requested threshold and possibly, there are also strong theoretical 
indications of existence of stabilization, meaning by this that the ionization 
rate goes to zero by increasing the intensity of the laser field 0. There 
has been a lot of theoretical work following the appearance of those effects 
and, as a general approach, people aims to develop models that account for 
the physics of strong-laser atom interaction trying to support them both by 
numerical and experimental work. So, a very succesful model for harmonic 
generation has been firstly put forward in [Q]. This is the recollision model 
where is assumed that the outer electron goes into continuum due to the 
laser field and here is accelerated by the field, this step being described classi- 
cally, then it recombines with the core generating the harmonics. A quantum 
improved version, understanding the emitted radiation as bremsstrahlung ra- 
diation due to the electron approaching the core has been given in This 
improvement has the advantage that an explanation of the appearance solely 
of odd harmonics of the laser frequency is given. Another quantum account 
of the recollision model has been given in Q. To obtain these quantum mod- 
els a number of reasonable assumptions have been made as the disregarding 
of intermediate states, the absence of resonances and so on, that need to be 
justified. These quantum versions of the recollision model prove to be a satis- 
factory explanation in all the experiments carried out so far with harmonics. 
Then, any theory starting from quantum electrodynamics has to cope with 
such a succesful understanding of the situation at hand. 

Since the initial work by Kulander et al. a lot of numerical work is 
currently performed to understand strong laser-atom interaction. Anyhow, 
although this work is highly interesting our aim in this paper is quite dif- 
ferent. That is, we are trying to rederive the recollision model from purely 
quantum arguments without any other approximation than the schemes of 
current experiments in high-intensity laser field require. Our conclusions con- 
cern the justification needed to understand the theoretical models discussed 
in Ref. [^,^. In this light, we take for granted that numerical computations 
support the recollision model as the experiments do. 

After the pioneering work in Q, a lot of studies have also been carried 
out assuming a simple two-level model. So, one of the open problems in this 
field of research is if such a simple model really accounts for the physics of 
harmonic generation. In this light, the general structure of the spectrum has 
been obtained by Floquet theory and a parallel with the recollision model 
has been given in [0,^ for the so called cut-off law, that is the rule to obtain 
the region of the spectrum where the intensity of the harmonics goes rapidly 
down. It should be said that such a model cannot account for ionization and 
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so appears somewhat rough at best. 

In the hght of the above approaches, in this paper we will try to change 
the point of view beginning directly from quantum electrodynamics, taking 
as a starting point our work in Ref. Q that here is deepened and improved. 
Indeed, modelling is an important way of understanding but, if one is able to 
solve the equations of the full theory without resorting to numerical methods, 
again a deep way to understand physics is also given. So, our main aim is 
to study in depth what is the physics of strong laser-atom interaction by 
directly solving the Schrodinger equation. The theory we obtain is anyway 
enough general to be possibly applicable to other physical situations where 
strong electromagnetic fields play a role as e.g for Rydberg atoms. 

It is interesting to note that, in view of using the dual Dyson series pre- 
sented in [10 1 and discussed here by using the Birkhoff theorem [12|, the Pauli- 
Fierz transformation and its classical counterpart, the Kramers-Hennerberger 
transformation, turn out to play a dominant role. Indeed, when unitary 
transformations are applied, using probability amplitudes permits to extract 
physics without worrying about. This is a standard approach e.g. in quan- 
tum field theory for the interaction picture. Things do not change for the 
dual Dyson series through the Pauli-Fierz transformation. 

Some approximations to start with are needed, but we take them directly 
from the experiments to tailor the Hamiltonian of quantum electrodynamics 
to the kind of problems we want to discuss. So, we take the electromagnetic 
field to be second quantized and the atom field in the non-relativistic approx- 
imation. The quadratic term of the vector potential in the Hamiltonian is 
considered just to see, through a single mode approximation, that taking it 
into account would shift the frequency of the laser in the harmonics. So, the 
condition to neglect it is also given. Finally, as the wavelengths of interest are 
much larger than the atomic radius, the dipole or long wavelength approxi- 
mation is also taken. Then, a Pauli-Fierz transformation is applied to study 
strong laser-atom interaction. We will show that only if a large number of 
photons is present and the amplitude of the quiver motion of the free-electron 
is much larger than the atomic radius, a condition normally met in this kind 
of experiments and a key approximation that we will assume in this paper, 
than just odd harmonics are generated, assuming the atomic potential to be 
spherically symmetric. Indeed, for a large number of photons the Pauli-Fierz 
transformation reduces to the Kramers-Henneberger transformation [^] as 
should be. 

An interesting result obtained in this way is that, the atom appears to 
be a kicked quantum system, being this result rigorously derived from the 
Hamiltonian of quantum electrodynamics. The kicking hypothesis has been 
assumed for Rydberg atoms in a microwave field firstly in Ref. ||l^ in order to 
explain the experimental results about ionization. Indeed, kicking can mean 
localization as a counterpart of classical chaos [16|. This is one of the main 
results of this paper pointing toward the merging of such different fields of 
research. 
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A kicked Hamiltonian can in principle be solved exactly. Anyhow, we want 
to study an atom in a strong laser field doing perturbation theory. So, we 
will obtain a perturbation series having the product u;7 between the Keldysh 
parameter 7 = ^ being Ib the ionization energy and Up = the 
ponderomotive energy for an electric field E, and the laser frequency cv smaller 
than the distances, properly shifted by the laser field, between the energy 
levels. This accomplish the task to obtain a perturbation theory dual to the 
standard time-dependent theory applied to an atom in a weak electromagnetic 
field. Using this approach, we are able to show that the rate of above threshold 
ionization determines the duration of the harmonics in the spectrum and then, 
also the form of the spectrum itself by Fourier transform. Beside, if a rough 
model for the outer electron in the field of the rest of the atom is taken having a 
Coulomb form with a Z^ff for the atomic number to give the correct ionization 
energy, it can be seen that the duration of harmonics agrees fairly well with 
experimental results as given e.g. in |17]. This simply means to scale all 
the formulas for hydrogen-like atoms by ionization potential, ponderomotive 
energy and laser frequency. 

The applicability of perturbation theory is indeed possible as, through 
a modified Rayleigh-Schrodinger perturbation theory one can show that the 
energy levels of the atom are shifted in such a way to change very little the 
wave function giving a prove of stabilization in the limit of laser frequency 
increasingly large. This is what we call "rigidity" of the wave function, very 
similar to the behavior of a superconductor due to the presence of the binding 
energy of the Cooper pair. This should give a hint, in the framework of 
quantum chaos, to study the change of statistics of such energy levels, possibly 
by some numerical work. 

Let us finally point out that a number of methods have been devised 
starting from Volkov states, that are the solution of the Schrodinger equation 
for a free particle in a plane wave. The prototype of these approaches is 
given in Ref. |18|. Anyhow, our point of view assumes different asymptotic 
states that are naturally derived from dual Dyson series as an entangled state 
between the radiation and the atomic state: This is just the leading order 
approximation. Then, e.g. no hypothesys a priori on the atomic potential 
being short ranged is needed. On the other side, no general proof is known of 
Volkov states being the proper asymptotic states for perturbation theory in 
a strong electromagnetic field. Beside, it appears very difficult to understand 
why just odd harmonics are experimentally observed as this approach keeps 
both odd and even harmonics Anyhow, an improved version has been 
recently given in Ref. by the Volkov states in a second quantized field 
giving just odd harmonics. This approach satisfy the principle of duality in 
perturbation theory but what is changed is the initial state used to obtain the 
perturbation series and is generally presented as a non-perturbative method 
through the results for the computation of the amplitudes given in Ref. |18|. 

The paper is so structured. In sec.|| we give a brief presentation of dual- 
ity for the Schrodinger equation and discuss the Birkhoff theorem using the 
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dual Dyson series. In sec. Ill we give the Hamiltonian formulation of quantum 
electrodynamics to start with and show why the quadratic term is negligible. 
Then we perform a Pauli-Fierz transformation to obtain the dual Hamiltonian 
and show how it reduces to the Kramers-Henneberger form. A discussion of 
the two- level approximation with respect to the full theory is also given. In 
sec. IV we apply the approximation of large amplitude of the quiver motion 
in the laser field with the respect to the atomic radius to show how kicking 
arises and why the atomic wave function turns out to be rigid in the sense 
given above. In sec.^ perturbation theory is done in the tunnelling regime de- 
riving the rate of ionization and consequently the spectrum of the harmonics. 
Finally, in sec. VI the conclusions are given. 



II. DUALITY AND BIRKHOFF THEOREM 

The duality principle in perturbation theory, when applied to the 
Schrodinger equation written as (here and in the following h = c = 1) 

states that, by choosing as unperturbed Hamiltonian, the perturbation 
series has a development parameter exactly the inverse of the series obtained 
by taking Hi as unperturbed Hamiltonian. Formally, this latter case means 
that we are multiplying Hi by an ordering parameter A going to infinity. 



Duality principle is enough, as shown in Ref. [10,1^, to prove that the dual 
to the Dyson series is just the adiabatic approximation and its higher order 
corrections. Indeed, from the theory of singular perturbation we can rederive 
a similar result for the equations of probability amplitudes by the Birkhoff 



theorem [12|, where "singular" means that the small parameter multiplies the 
higher derivative in the equation. In this way, a dual Dyson series can be 
obtained also for the probability amplitudes. In fact, the dual perturbation 
series is obtained by solving at the leading order the equation 

Ai/i(i)|V^(o)) = i 



,(0)v_.5|V'(°)) 



dt 

taking A — > oo. But this can be stated as a singular perturbation problem. 
Indeed, using the eigenstates of the unperturbed part HQ\n >= En\'n > and 
putting \ip{t) >= Y^n^~^^"^C"n.{t)\n >, one has for eq.([l|) the system of linear 
differential equations for the amplitudes 

.^dorn^ = ^e-(^"-^-)*(m|i7i(t)|n)a„(t) 

n 

being now ^ = j- If the Hamiltonian Hq has a finite set of N eigenstates and 
eigenvalues and considering the amplitudes am{t) as elements of a vector a(t), 
we can apply a result obtained by Birkhoff on 1908 in singular perturbation 



theory |12], that states that the equation (|3D has a fundamental set of sohitions 
hj{t) with 1 < J < if the eigenvalues Xj{t) are distinct, given by 

rt 



bkj{t) = exp 



-I 



to 



Xj{t')dt' 



ex-p[i'jj{t)]ukj{t) + 0(e) 



k-th element of the vector hj{t), Jj{t) = f^^Uj{t')i-^\ij{t') the geometrical 
part of the phase, Uj{t) the eigenstates of the matrix h aving elements A^nn — 
e~'^^^"~^"^^^ {m\Hi{t)\n) . The condition of no crossing of eigenvalues Xj{t) is 
also required. The final leading order approximation is then written as 



N 



i=i 

with the coefficient aj given by the initial conditions a(0). Again, we realizes 
that formally this is the adiabatic approximation applied to the operator 



e'-^o*i?i(t)e-*-f^o*. Then, we are arrived at similar conclusions as in Ref. 
for dressed states where, instead, duality principle was used, in the limit 
A ^ oo. But, the use of duality principle permits to get rid of the limitation 
on the spectrum of Hq being bounded to eigenstates and eigenvalues and 
then, it appears as a more general tool to treat also this kind of problems. 
So, we relax the condition on the spectrum of Hq required by the Birkhoff 
theorem, as duality permits to give an alternative way to derive it. 

Indeed, one can build the dual Dyson series by also doing an unitary 
transformation that removes the perturbation Hi{t) in eq.([l|). This defines a 
dual interaction picture. In fact, one can see that the unitary transformation 
operator must be a solution of 

and so Upit) turns out to be the evolution operator of the adiabatic approx- 
imation in the limit A — > oo, proving the equivalence. But, this equation 
can also be solved exactly in some cases, without resorting to the adiabatic 
approximation giving the series 



IV'(t)) = f/f(t)rexp 
with T the time-ordering operator 



mo)) 

to 



III. NON-RELATIVISTIC FORMULATION OF QUANTUM 
ELECTRODYNAMICS 

A. Full Hamiltonian and the Quadratic Term 

The starting assumptions about the formulation of quantum electrody- 
namics that we need are the following. Firstly, we takes a non-relativistic 
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approximation and neglect spin effects. Secondly, we cut-off the wavelengths 
limiting the analysis to the long ones that is, we take the dipole approxima- 
tion: currently, this is in agreement with the experimental results for harmonic 
generation. Finally, we do second quantization on the electromagnetic field 
as we want to understand the properties of the scattered light in experiments 



with strong fields. So, we write the Hamiltonian in the Coulomb gauge |13| 
(here and in the following h = c = 1) 



2 



H = Y^ a;kai(k)aA(k) + |- + ^(x) + - Ap + ^ 

2m m 2m 



A,k 



with A = 1,2 meaning a sum over polarizations, V^(x) the atomic potential, 
a^(k),aA(k) the creation and annihilation operators for the mode k with 
polarization A. So, we have 

A,k 

being e;^ the complex polarization vector, = \fy^ ™^ ^ normalization in 
a box of volume V is assumed everywhere. 

To see how much relevant is the quadratic term, we refer to the experi- 
mental result that just harmonics of the laser frequency are observed. Then, 
specializing the above Hamiltonian to a single mode we derive the Heisen- 
berg motion equation for the creation and annihilation operators. So, we take 

2 2 

H = ua^a + — + y(x) + -Ap + —A\ 
2m m 2m 

Assuming a linear polarization, that is 



— e(at + a) 

we can use the same argument given in Ref. |21] that the laser frequency 



should be shifted by squeezing. Indeed, one has the Heisenberg equations for 
a and 

da e^g^ 4. eg 

— = —iLoa — I (a + a') — I — ep 

at m m 

da'^ e^g"^ x eg 

— — = luja + I (a + a^) + i — ep 

dt m m 

and without going into details of computation, by duality, as a first approx- 
imation, one neglects the atomic Hamiltonian and then, managing p as a 
c-number, the time dependence of the creation and annihilation operators is 
harmonic with frequency J7 = \J uj"^ + and not just u! as should be required 
by the experimental results. So, squeezing terms can be neglected otherwise a 



shifted frequency of the harmonics would be observed |21]. This effect should 
be more pronounced as greater is the order of the harmonic. But, consid- 
ering this shift with current frequencies and gas densities really negligible. 



the quadratic term of the field in the Hamiltonian (^) can be systematically 
neglected giving the final Hamiltonian 

2 

H = J2u;^al{k)ax{k) + ^ + ^(x) + ^ ^^^(e^Kk) + exax{k))p. (13) 
This will be the starting point for our further analysis. 



B. Pauli-Fierz Transformation and Kramers-Henneberger Hamiltonian 



We now try to approach the Hamiltonian ( p^ ) using duality principle, 
that is, we try to compute the dual Dyson series. So, one can remove the 
perturbation by a Pauli-Fierz transformation given by [l^] 



UpF = exp 



5](/3I(k)a;,(k)-/3;,(k)4(k) 

k,A 



where one has 



(14) 



/?A(k) 



-5ke -P 



and obtains the Hamiltonian |13| 



HPF = E ^k4(k)aA(k) + ^^+V 

k,A 



x-i5](/3Kk)aA(k)-/3,(k)4(k 



k,A 



being m* the renormalized mass due to the field. Here we try another way, 
to agree with duality in perturbation theory. That is, firstly we consider the 
Hamiltonian (|l3|) in the interaction picture giving 



2 

Hi = ^ + y(x) + - 5] 5k [e^4(k)e^"''* + eAaA(k)e- 

A,k 



P- 



Then, we solve the leading order equation 



e 

m 



A.k 



The solution is standard and can be written as 



Upit) = exp 



E7k,A(t)(ei-p)4(k) 

k,A 



exp 



E7k,AW(eA-p)aA(k) 
k,A 



exp 



(15) 



(16) 



(17) 



(18) 



with 



e / 27r e^'^kt _ i 



7k,A = — 



E"k,A(i)(eI •p)(eA -P) 
k,A 

(19) 



(20) 
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and 

Ok.A = -^l7k,A(t)P- 

So, we have found in this way a time-dependent version of the above Pauh- 
Fierz transformation. But now, by eq.(0) we have a dual Dyson series for 
quantum electrodynamics. From this, one has e.g. that the leading order 
wave-function in the interaction picture, if we assume the atom in the ground 
state, is given by 



\i;('\t))j = Upms)\ 



a) 



being the atomic state and \a) the initial state of the laser field described 
by a coherent state. As it should be expected this is an entangled state 
between the atomic state and the field state It is easy to realize, using 
the momentum representation for the atomic state, that the scattered light is 
indeed coherent and preserves the property of the initial state for the field. 
On this basis, we expect the light of the harmonics to be coherent as well 



21 1 . In fact, in order to evaluate higher orders in the dual Dyson series we 



rewrite the Pauli-Fierz transformation as 

Upit) = exp J2 (7k,A(*)(eA • p)aA(k) - 7k,A(i)(eA • p)4' 

k,A 

obtaining the transformed Hamiltonian for higher order computations 

Hf = ulit) + y(x)^ Upit) = 1^ + y [x + x(t)] 

being 

Mt) = -^E (7k,AWeAaA(k) -7k,A(t)6l4(k) 



k,A 

Now, we specialize the above construction assuming that initially the atom 
is exposed to an intense laser field. It is well-known that the creation and 
annihilation operators can be expressed through the number operator rik^A = 
A^k,A and a phase operator by the Susskind and Glogower or the Pegg 
and Barnett construction |2^. But, the distribution probability of the phase 
can be made to coincide in such a way to give for a coherent state a definite 
phase having quantum fluctuactions that go to zero for a very large mean 
number of photon, that is the case of the intensity of the laser field of interest. 



Then, one can keep just one mode and operate the substitution |25] 



The Hamiltonian Hp is then reduced to the Kramers-Henneberger form |14| 
for a classical field that is, 

Hhk = ^ + V[x- XLe{sm{iot -</>)+ sin(</.))] 
2m 
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where a linear polarization has been assumed, all subscripts have been 
dropped and A^, the free-electron maximum excursion in a monochromatic 
field, can be rewritten as Ai, = being E the intensity of the laser field. 
The above Hamiltonian can be obtained by an unitary transformation, the 
Kramers-Henneberger transformation, for an atom in a monochromatic clas- 
sical field. Then, we simply neglect quantum fluctuaction beyond the leading 
order due to the characteristic of the electromagnetic field. The relevance of 
this limit for a coherent state is that it makes consistent the arguments that 
follow about the properties of the harmonics. We just note that the advantage 
of such a derivation is that it gives us the coherence property of the scattered 
laser light in atomic scattering experiments, and permits to understand the 
way from a quantum electrodynamics formulation to a classical Hamiltonian 
that fully accounts for the situation at hand. 

As a general approach, in agreement with the duality principle in pertur- 
bation theory , we see that to study strong laser-atom interaction we apply 
unitary transformations to remove the big part of the Hamiltonian that in this 
case is due to the field. We now apply the same idea also to the two- level 
model going to a dual interaction picture. 



C. Two-level Approximation 

The two- level model generally adopted in the study of harmonic generation 
can be cast in the form (e.g. |^) 

A 

H = "^^3 + cos{Lot) 

being A the distance between the two-level, Q the strength of the field, cv the 
frequency of the laser field and ai, the Pauli matrices. The first point to 
note is that, in order to be in agreement with the recollision model, one of the 
two levels should be in the continuum or, at best, ^ is the ionization energy Ib 
0. In this way, one can rederive the well-known cut-off law for the frequency 
of the maximum harmonic Ib + 3C/p, being Up = IrnL'^ ponderomotive 
energy that in this model should be proportional to Vl as pointed out in 
Ref. 0. But, in this model the continuum is missing and then, ionization 
cannot be described. This means to lose any connection between ionization 
and harmonic generation that, as we will see, plays a role in the duration of 
the harmonics itself. Finally, let us apply a Pauli-Fierz transformation (or 
Kramers-Henneberger transformation) by removing the oscillating term as 



Upit) = exp 



-zo"i— smfwt) 

UJ 



to obtain the transformed Hamiltonian 
J Jo ( — j + Acj3 J2n [ — j cos(2nwt) + Aa2 J2n+i [ — j sin((2n + l)Lot). 

n=l n=0 
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From this Hamiltonian is quite easy to realize that both odd and even har- 
monics play the same role and the only way to make ones or the other appear 
is to choose proper initial conditions. To show that things are really in this 
way, we turn back to Hamiltonian ( p8|) and apply the Birkhoff theorem to 



rederive the spectrum already obtained in Ref. |11|, given in a general form 
by the Floquet theory in Ref. P . 
So, we have to solve the equation 



A 



(T3 + OO"! COs(wt) 



\m) 



: d\m) 

dt 



that in the interaction picture becomes the system of equations 



.d_ / ai(t) 
dt I 02 (t) 



ne-'^^ cos{ujt) \ ( ai{t) 

fte'^^ cos{Lot) ] I 02 (t) 



to which the Birkhoff theorem can be applied. It should be pointed out 
that this approximation is consistent with A ^ and we take uj <^ A, 
in agreement with the identifications of Ref. and in the way the model 
should compare with experiments. Otherwise, the Birkhoff theorem should 
be applied differently. Under these conditions, we show that even harmonics 
can be present depending on the initial conditions. These harmonics are 
known as hyper- Raman lines in the current literature |^,26|. 

Firstly, we note that eq.(^) is the same as eq.(pTD when use has been 
made of the eigenstates of a^. Then, we can simply use the results of Ref. |11| 
about the dressed states of Hamiltonian (Bsl) and obtain the Birkhoff basis as 



and 



bl(t) 



b2(i) 



e -2-6^' 



V2 



'iAt 



^/2l 1 



6 2 6 



Then, we conclude that the Birkhoff theorem gives the same results obtained 
by other approaches in Ref. [ p^j2^ . The spectrum is then given by [11 1 



m)Wi\m) « a2(0)at(0)e-*^«* + a^(0)ai(0)e*^«* 



a2(0)P)AE^....i(g^) + 



n=0 



+ ^(a^(0)ai(0)e^^«* - a2{0)aU0)e~'^^')A ^ J2„ f — 



n=l 



2i7\ s'm(2nujt) 



with A/j = AJo [j^^ the renormalized level separation. This form of the 
spectrum is in agreement with the one derived by the Floquet theory as given 
In the range of validity of this approximation as stated above, this 



m 



proves the assertion that the even harmonics should be considered on the 



same ground as the odd ones for the two-level model. As we are going to 
show, this makes the applicability of this approximation to harmonic genera- 
tion somewhat unappropriate as the appearance of just odd harmonics in the 
experiments has a deep physical meaning that here is overlooked. In fact, the 
two-level model relies just on the initial conditions to select the properties of 
the spectrum. 



IV. PROPERTIES OF ATOMS IN INTENSE LASER FIELDS 

A. Odd Harmonics and Kicking for an Electron in a Coulomb Field 

We now turn our attention to the classical limit of the Pauli-Fierz Hamil- 
tonian that is, the Kramers-Henneberger Hamiltonian that here we rewrite 
assuming the phase (j) = obtaining 

p2 

Hhk = 7, h 1^ [x - Alc sm(tjt)] . 

The potential has the Fourier series ^ 

oo 

V[x- XL€sm{ujt)] = J2 Vk{x)i'' k'^^* + (-l)"e-*'^'^* 

k=0 



where 



Vk{x)= V{x-eXLx) 

J-i vr V 1 - x'^ 



and being Tk{x) the A;-th Chebyshev polynomial of first kind. In this sec- 
tion we study the terms having k > 1, deserving the analisys of the static 
contribution with fc = 0, i.e. the dressed potential, for the next section. 

We now take the unperturbed potential l^(x) to be spherical symmetric. 
It easily realized that one has different parity for even components being 
V2n{-^) = t^2n(x), and odd components being 't;2n+i(-x) = -V2„+i(x). But, 
the unperturbed part of the Hamiltonian is invariant by parity transformation 
while, the odd time-dependent part breaks this simmetry. So, by perturbation 
theory one has at first order {iij{t)\x\ij{t)) = (V'^^HOklV'^^H*)) + cc. being 
I'i/^^'^)) the initial atomic state such that {il)^^\t)\x\il)^^\t)) = 0, and the only 
way not to obtain a null spectrum through perturbation theory is from the 
odd components. This result turns out to be in agreement with the general 



one obtained through Floquet states in Ref. |27]. 

To apply the above scheme to the Coulomb potential one has to consider 
the integrals 

, , Ze2 /-i dx' 1 Tkix' 



+KXZJ 

but we have to study them around the origin of coordinates where the 
Coulomb potential turns out to be singular. Indeed, at small distances the 
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theory given above has its shortcomings due to the dipole approximation and 
the neglecting of the relativity. Then, some regularization is needed for the 
above integrals. So, let us study the integrals (^) at the point y = z = ao, 
introducing in this way the cut-off ?? = <C 1, and we are reduced to the 
analysis of the leading term 

Ze^ fi dx' 1 Tk{x') 



x' + 27]^ 



very near the origin. These integrals turn out to be zero for odd Chebyshev 
polynomials when one takes x = and, for even Chebyshev polynomials one 
gets an extremum such that the integral is bounded. The extremum depends 
on 7] but increases slowly. Then, we can remove the even terms, after we take 
V2n at their values at the origin, by an unitary transformation. The same 
cannot be done for odd terms as they go to zero at the origin. These integrals 
can be studied by a Taylor series at x = giving for the dipolar term 



dip 



X 



Ze^ dx' 1 Tk{x')x' X 



Al J -I vr (x'2 + 2r/2)f Vl - x'2 Al 
and the integral gives the cut-off dependent approximation 



1 dx' 1 Tk{x')x' 2 , , ^. ^r, 

— T7> ^ A^Ao ^ — log r/ 2n + 1-1-. 

1 ^ {x'^ + 2rf )2 Vl - x'^ vr 



that has a logarithmic divergence for the cut-off r/ going to zero. If we intr- 
douce the ratio ^ ^ 1 that defines the region where we want to study the 
physics of the model, we can reabsorb the divergence — log(r/) as a renormal- 
ization into it and introduce the physical ratio ^ = This means that we 

have two bare constants eg and rj and this goes to zero through the relation 
= e "0 with e the physical ratio. Then, the limit eg ^ implies r/ — > 0. 
Through this redefinition, the dipolar terms are 

^f(x) = -6^(-ir(2n + l)^ 

In this way, we can now prove that the electron in a Coulomb field and an 
intense laser field, such to have the parameter e ^ 1, undergoes kicking 



turning the problem into one of quantum chaos |16|. 

The question of a kicked hydrogen atom has been put forward to explain 
experiments on ionization of Rydberg atoms in an intense microwave field 
through the mechanism of dynamical localization If]. Pioneering work 
has also been done on strong fields and atoms, using kicked models e.g. in 
Ref. But, a rigorous proof of this behavior derived directly from quantum 
electrodynamics has not been given yet. Here, we accomplish this task by com- 
pleting the analysis of the time-dependent part of the Kramers-Henneberger 
Fourier components. To complete this derivation, we point out that a periodic 
distribution as ^^(t) = J2t=-Qo^(.'^ ~ that has period T, can be defined 
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through the coefficients of its Fourier series. So, e.g. Sxit) = ^ J2t=-oD e*^"'^T 
has all constant coefficients given by To have convergence of the Fourier 
series in the sense of distributions, the coefficients fn must satisfy the cri- 
terium of being slowly varying [30|, i.e. \fn\ < Mlnl*^, being M and k two 



constants. 

Then, by using eqs.(p3|) and one gets the dipolar potential as 



,H = -i=^—e V (2fc + l)e^(2fc+iM ^2e4^— V(2/t + l)sin((2fc + l)a;t). 

«B k=-oo as 

The Fourier series appearing here has convergence just in the sense of dis- 
tributions. We now give an explicit form for the dipolar term. Indeed, one 
has 

+ 00 y rp\ r) +00 



2 



k=—oo 



2 T 



so that, the dipolar term is defined through the derivative of a periodic dis- 
tribution of period ^. One can write at last 

TT Ze^ X d ~ , . 

VKH = -e—- ^ST{t) + --- 

uj'^ \l as dt 2 

and we draw the conclusion that the motion of an atom in an intense laser 
field undergoes kicking with period ^ so that, localization, a typical effect of 
quantum chaos, can happen. 

Now, we show that the Rydberg atoms used in ionization experiments 
with microwave are indeed in a regime where the ratio between the atomic 
dimensions and the amplitude of the quiver motion of the free-electron in 
the laser field is largely smaller than unity, and the above theory applies. As a 
starting point we take Ref. for a typical experiment. The less favourable 

case is for an intensity of the microwave field of 2.5 V/cm, at 12.4 GHz and 

2 2 

the Rydberg atom has nn=98. This gives for the ratio "9"^ = the 

value 0.0027, largely lesser than unity. Instead, the most favourable case is 
given by the intensity of the electric field of 21 V/cm, at 18 Ghz and no=64 
giving a ratio of 0.00029 improving the situation of about a magnitude order. 
Then, we can conclude that Rydberg atoms are kicked by the microwave field 
and localization can happen, in agreement with all the current researches. 
This result comes directly from quantum electrodynamics and so, it is fully 
justified. 

For the case of the harmonic generation the situation is still better because 
the Keldysh parameter 7 given by \J -^ff- can be lesser than one permitting 
a straightforward application of the perturbation theory. A typical example 
of an experiment in this regime is given in Ref. This possibility can be 
fully exploited if a set of unperturbed states can be found to do perturbation 
theory. This is indeed the case as we will see in the next section. 
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B. Rigidity of the Atomic Wave Function 



In order to do perturbation theory, one generally needs a full set of or- 
thonormal functions to start with, representing the unperturbed system. But, 
here one has to diagonalize the Hamiltonian 



2m J-i IT y/T 



X' 

that is generally an impossible task analitically. 

What we want to do is to analyze this problem in the limit where the 
amplitude of the quiver motion of the free-electron in the laser field is much 
larger than the atomic radius, further specializing the above problem to the 
Coulomb potential. 

To evaluate the degree of deformation due to the laser field in this approx- 
imation we apply a modified Rayleigh-Schrodinger approximation as obtain- 
able from the time-dependent perturbation series. This can be accomplished 
easily if one makes the multipolar expansion of the dressed Coulomb potential 
in eq.(^) as 

2n 



n=l ^ ' ^ ' 



being An = j\ dxx^^ / {iry/l — x^) and P„ the n-th Legendre polynomial. 
So, we approach the static part of the Kramers-Henneberg Hamiltonian dif- 
ferently from the time-dependent part. The reason to do that is that we 
expect very large shifts of the energy levels and that just very few terms of 
the multipolar series ( ^8| ) really contributes to the matrix elements in the 
Rayleigh-Schrodinger series. 

The way we compute the Rayleigh-Schrodinger corrections is taken from 
the time-dependent perturbation theory. Indeed, let us suppose that one can 
neglect the time dependent part of the Kramers-Henneberger Hamiltonian, 
so that one is left with the time-independent problem 



where we have put 



Ha = ^- — + 5Vkh{^) 
Zm r 



n=l 



r J ^" \ r 



and we want to solve the problem 

This is the way stabilization is studied in the Kramers-Henneberger frame 
when the limit u; — > oo is taken [^,32|. Here we are able to prove that stabi- 



lization indeed exists in the limit of large ratio of free-electron quiver motion 
and Bohr radius due to the rigidity of the atomic wave function. 
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We take as unperturbed states the ones of the Coulomb problem by setting 

m)) =T.^n{t)e-^^An) 

n 

being 

One gets for the amplitudes 
%a^{t) = {m\SVKH{y^)\m)amit) + ^ e-^(-^"-^-)*(m|,5FKi/(x)|n)a„(t) 

and introducing bm(t) = e~"^^'"*am(t), SEm = {m\SVKH{'x.)\m) and Em = 
Em + SEm, we arrive finally at the equations 

ihm{t) = J2 e-'^^"-^-^'{m\6VKH{^)\n)bnit). 

Then, if the shifts SEm are really large in the limit of large ratio between the 
quiver amplitude of the motion of the free electron and the Bohr radius, the 

Rayleigh-Schrodinger corrections to the initial wave-function are really small 
and this state is "rigid" with respect to the perturbation introduced by the 
laser field. Indeed, one has 

bmit) = 6.(0) + E e--(^"-^-)* ^"^'f"_"|^^'"^ 6„(0) + . . . 

where an adiabatic switching of the perturbation has been introduced. We see 

that the first order correction gives a modification of the Rayleigh-Schrodinger 
perturbation theory as in place of the unperturbed energy levels Em there are 
the modified energy levels Em- This appears also as an improvement with 
respect to the Brillouin-Wigner perturbation series. Indeed, one has for the 
wave function 

n n k^n En 

In this way is possible to show that the first-order correction is really small. 

For the level shift one gets for the first few levels (we have put as for the 
Bohr radius, I for the orbital angular moment and Iz for the third component 
of the orbital angular moment): 

n=l 

(5£;ooi = (/ = 0, = 0, n = l\5VKHi?^)\l = 0, = 0, n = 1) = 
n=2 
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(5^102 = (/ = 1,1, = 0,n = 1\bVKH{^)\l = 1,1, = 0,n = 2) = — ^ 

(5£;ii2 = bE^^vi = (Z = 1,Z, = l,n = 2|,5F;^,j(x)|Z = 1, = 1, n = 2) = — ^eV^L^ ' 



480 as VflB, 

<5£;oo2 = (Z = 0,Z, = 0,n = 2|5y^/^(x)|/ = 0,Z, = 0,n = 2) = (59) 



n=3 

(5^223 = (/ = 2, = 2, n = 3|<5Fk//(x)|Z = 2, Z, = 2, n = 3) = -— ( -- ( ^^^^ ) + 

Ob \ 57dU \aB / o1d48U \aB J 

dE2is = = {l = 2,l, = l,n = 3\SVkh{^)\1 = 2, Z, = 1, n = 3) = — ( f — V + ^ ^^"^ 



^eV 1 /Al^ 1 (XlV 



as yll340 VaB/ 204120 V as 

^^^^ '-(^W^f^Y] 



ttB \ 5670 Vas/ 136080 V«_B/ 
(5^103 = (Z = 1,Z, = 0,n = 3|,5FKff(x)|Z = 1,Z, = 0,n = 3) = -i-— (^Y 

(5^113 = 5^1-13 = (Z = 1, Z. = l,n = 3\6Vkh{^)\1 = l,h = l,n = 3) = -—- 

1620 qb \aB J 

5Eoos = (/ = 0, Z, = 0, n = 3\6Vkh{^)\1 = 0, Z, = 0, n = 3) = (60) 
n=4 

<5-E334 = (Z = 3,Z, = 3,n = 4|<5Vk//(x)|/ = 3,Z, = 3,n = 4) 

Ze^ f 1 ^Al\' ^ 3 /XlV ^ 25 /Al 



as y 32256 V OB 7 50462720 V as/ 113359454208 Vofi 
(5^324 = SE3-24 = {l = 3,l,=2,n = 4\SVkh{x)\1 = 3, Z, = 2, n = 4) 
Ze^ 1 ^A^y ^ 25 fXL 



as \^ 7208960 V OS/ 18893242368 Vob 
SEsu = SEs^u = {l = 3,h = l,n = ^SVKni^W = 3,Z, = l,n = 4) 

^ / l_ (XlV 1 125 /AL^6^ 

as V 53760 Us/ 50462720 U-B / 37786484736 Ub/ 
SE304 = {1 = 3,1, = 0,n = A\5Vkh{x)\1 = 3, Z, = 0, n = 4) 

Ze2 / 1 /Ai\2 3 /AlV 125 f 



as \^ 40320 VobJ 25231360 Vob 7 28339863552 Vofi 
5^224 = 5-E2_24 = (Z = 2, Z, = 2, n = 4\SVKHi^)\l = 2, Z, = 2, n = 4) 

aB \13440 Vafi/ 4587520 Vos/ / 
dE2u = SE2-U = {l = 2,h = l,n = 4|(5y^ij(x)|Z = 2,Z, = l,n = 4) 

as 1^26880 Ub/ 1146880 Us, 
5^204 = (Z = 2, Z, = 0, n = i\5VKHix)\l = 2, Z, = 0, n = 4) 

UB V 13440 Vasy 2293760 VaB/ / 
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8E 



114 



()^i_i4 = {l = l,h = l,n = 4|,5Vi^H(x)|/ = 1,/, = l,n = 4) 



1 Ze2 /A 



<^£^104 = (Z = 1,/, = 0,n = 4|^FKi/(x)|Z = 1,Z^ = 0,n = 2) 



1 



3840 Ob Ub. 
Ze2 /\-x2 



1920 as \aBj 



SEooA = {1 = 0,1, = 0,n = 4\6VKH{y^)\l = 0,1, = 0,n = 4) = 

It easily realized that just very few terms of the multipolar series of 5Vkh 
really contribute to level shifts making the argument working. Beside, as it 
should be expected the s-states have no shift. 

Indeed, for the matrix elements, assuming as initial state the ground state 
of the atom one obtains 

n=2 



{l- 


= l,lz 


= l,n = 2\SVKH{yi)\l-- 


= 0,1, 


= 0,n = 


-1)=0 


{l- 


= l,lz 


= 0,n = 2\SVKH{^)\l-- 


= 0,1, 


= 0,n = 


-1)=0 


{l- 


= 0,1, 


= 0,n = 2\SVKH{:><i)\l-- 


= 0,1, 


= 0,n = 


-1)=0 


=3 

= 2,1, 


= 2,n 


= 3\SVkh{:x.)\1 = 0,1, 


= 0,n 


= 1) = 


Ze^ 1 
as 720 


= 2,1, 


= l,n 


= 3\SVkh{-^)\1 = 0,1, 


= 0,n 


= 1) = 





= 2,1, 


= 0,n 


= 3\SVkh{^)\1 = 0,1, 


= 0,n 


= l) = 


Ze'^ 

as 10800 


= l,lz 


= l,n 


= 3\6Vkh{^)\1 = 0,1, 


= 0,n 


= 1) = 





= l,lz 


= 0,n 


= 3\6Vkh{^)\1 = 0,1, 


= 0,n 


= 1) = 





= 0,1, 


= 0,n 


= 3\6Vkh{^)\1 = 0,1, 


= 0,n 


= 1) = 






n=4 

= 3, 1, 
= 3, 1, 
= 3, 1, 
= 3,1, 



3,n = A\6VKH{^] 
2,n = A\5VKH{^] 
l,n = A\6VKH{^] 
0,n = A\6VKH{^] 



= 2,l, = 2,n = A\5VKH{^) 

= 2,l, = l,n = A\5VKH{y^) 

= 2,l, = 0,n = 4\5VKH{:>^) 

= l,l, = l,n = A\6VKH{^) 
= l,l, = 0,n = A\6VKH{^) 
= 0,l, = 0,n = A\5VKH{-^) 



0,1, 

0,lz 

0,/. 



0,n 
0,n 
0,n 
0,n 










0,/^ = 0,ra = 1) 

/ = 0, Z;^ = 0, n = 1) 

/ = 0,/^ = 0,n = 1) 

l = 0,l, = 0,n = l) 
/ = 0,/^ = 0,n = 1) 
/ = 0,^^ = 0,n = 1) 



Ze^ 13^150 [XlV 
as) 



QB 150000 



Ze2 
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A, 



^ .V 

UB 15000 KubJ 




0. 



Again we see that just very few terms of 5Vkh really contribute to the matrix 
elements. Beside, this contribution is smaller if not zero with respect to the 
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level shifts. Then, these results strongly support the statement that the first- 
order correction into the series ( |57|) is really small in the limit ^ S> 1 as it 
should be. This also supports stabilization when the limit of laser frequency 
going to infinity is taken, with the ponderomotive energy being constant as in 
this limit, the time-dependent part of the Kramers-Henneberger Hamiltonian 
can be neglected 0j3^. 

Finally, we conclude that a fairly good approximation for the unperturbed 
Hamiltonian for our aims is 

HA^Y.^n\n){n\ (65) 

n 

where just diagonal terms are kept and the off-diagonal terms 

H'a= J2Hi^\i^\^^KH\n) (66) 

are neglected. We take the Hamiltonian ( |65| ) for time-dependent computations 
in perturbation theory. 



V. PERTURBATION THEORY FOR ATOMS IN INTENSE LASER FIELDS 

The experiments carried out to produce harmonics have a small laser 
frequency with respect to ionization and ponderomotive energy, then also the 
time dependent components of the Kramers-Henneberger Hamiltonian need 
to be considered. 

Our aim is to make a perturbaion analysis of the Hamiltonian 

+ 00 

X 



H = ^En\n){n\ +2euo-i—^{2k + l)s\n{{2k + l)ujt) (67) 

n 

being 



fc=0 



the Keldysh parameter. So, the perturbation theory is applicable in this case 
just when 7 <C 1 that defines the so-called tunnelling regime for the electron 
in the laser field. But, the situation here is more favourable as we have just 
to require that the distance between energy levels is much larger than ujj and 
we are able to account for a larger number of experiments than the tunnelling 
regime would permit. 

We write the equations for the probability amplitudes as 



A:=0 



Then, a generic term, out of resonance, will be written as 
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r' 



am{t) = am(0) + 



zew7^a„(0) 

n 



m 



X 



as 

,-i{{2k + l)L0+En-Em)t _ 

{2k + 1)U + En - E„ 



k=0 
+ 



/ i{{2k+l)uj-E„+E^)t _ -I 

{2k + 1) = ^ + 

'\{2k + 1)U -En + Em 





X 












as 





This result shows that when there are two resonant states m and n in the 
laser field, then we have Rabi flopping with Rabi frequency given by 

^ = e7(2fc + l)a; 

determined by the {2k + l)-th harmonic of the laser frequency. This implies a 
significant modification of the spectrum of the harmonics with respect to the 
observed patterns. Otherwise, out of resonance amplitudes are really small 
due to the large shifts in the limit of large ^ ratio. So, the interesting case 
is the resonance with the set of continuous states of the atom representing 
the situation of the recollision model. In this case, instead of fiopping, we can 
have decay, that means ionization, and a probability for the electron to turn 
back to the core emitting radiation. 

Then, we can apply the Wigner-Weisskopf argument to eq. (|70D . The 
first approximation is to reduce the system to the two levels really resonating 
and to assume all the other amplitudes to be small in the sense given above. 
This means to approximate eq.(^) by 



do(t) 
hp{t) 



p 



-i(Ep-Eo)t 



CLB 



ea;7e"*(^o-i?p)t / p 



as 



p)ap(t)^(2A; + l)e*(2'=+iM 

i{2k+l)ujt 



k=0 



o)ao(OE(2^ + l)e 

k=0 



where we have labelled the eigenstates of continuum for the atom by the 
momentum p and Eq = —Ib- We seek a solution of this equations by setting 



ao{t) 



e 2 



By substituting eq.([73D into eq.([7^) one arrives at the condition 



r 



• 2 2 2 \ " 





X 










as 





2 +00 



N2l 



i[Ep-Eo-(2k+l)u}]t+i^t 



y{2k + i) 

to {2k + l)u; - Ep + Eq 



For small V this indeed reduces to a time- independent expression [34| 

2 +00 



— =ie-fuj 



2„,2, ,2 ^ 





X 










CLB 





/c=0 



V 



1 



(70) 



(71) 



(72) 



(73) 



(74) 



{2k + l)uj - Ep + E^ 







iTT5{Ep -Eq- {2k + l)a;) 

(75) 



with V meaning the principal value. So, finally, one gets the a.c. Stark shift 
of the ground state of the atom given by 
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with atomic number ^e// can be adopted for complex atoms in a strong laser 
field and the above result used also in this case. As shown in Ref. the 
order of magnitude is correct. The same is true for T given by 



256 u;2 



■7 



n=nQ 



(2n + l)uj 



Ib 



{2n + l)uj 



that for the experiment described in Ref. [|17| yields for helium .012 eV and 
for neon .01 eV improving the results given in Ref. This in turn means 
that the mean lifetime of the harmonics ranges from 52 to 68 fs and, by 
e-folding, after a time of few 10^ fs the harmonics disappear in agreement 
with the experimental results. The result is independent from the order of 
the harmonics. We expect that the value of T should be in better agreement 
with the experimental results when Up increases or being the same, when the 
Keldysh parameter 7 decreases marking the so called tunnelling regime. 



VI. CONCLUSIONS 

We have given a full theory of harmonic generation starting from non- 
relativistic quantum electrodynamics in the dipole approximation. A relation 
with quantum chaos is obtained through a kicked model that originates from 
the model at small distances where a regularization procedure is needed. Per- 
turbation theory can be done in the tunnelling regime, where the Keldysh 
parameter is small, as the atomic wavefunction turns out to be rigid due to 
the way the energy levels are shifted by the laser field. In this way a deep 
connection between ionization rate and harmonic spectrum is given. Finally, 
a fully understanding of the recollision model is obtained and this was our 
main aim. 
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